Detection of Galaxy Cluster Motions with the Kinematic Sunyaev-ZePdovich Effect 
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Using high-resolution microwave sky maps made by the Atacama Cosmology Telescope, we for 
the first time detect motions of galaxy clusters and groups via microwave background temperature 
distortions due to the kinematic Sunyaev-Zel 3 dovich effect. Galaxy clusters are identified by their 
constituent luminous galaxies observed by the Baryon Oscillation Spectroscopic Survey, part of the 
Sloan Digital Sky Survey III. The mean pairwise momentum of clusters is measured at a statistical 
significance of 3.8 sigma, and the signal is consistent with the growth of cosmic structure in the 
standard model of cosmology. 

PACS numbers: 98.52.Eh, 98.62.Py, 98.70.Vc, 98.80.Es 


Introduction. The growth of cosmic structure over the ' history of the universe inevitably results not only in the 
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formation of dense objects, but also in motions of these 
objects. Measurements of these motions have the poten- 
tial to provide both a valuable consistency check on the 
standard cosmological model, and also an independent 
route to constraining cosmological parameters and the 
nature of dark energy. 

In 1972, Sunyaev and Zel’dovich realized that a moving 
galaxy cluster, which is largely composed of hot, ionized 
gas in a dark matter potential well, will induce a small 
brightness temperature shift in the microwave radiation 
passing through it. The shift is proportional to both 
the mass in electrons and the line-of-sight velocity of the 
cluster with respect to the microwave background rest 
frame [1, 2]. This kinematic Sunyaev-Zerdovich (kSZ) 
effect is distinct from the thermal SZ (tSZ) effect, in 
which scattering from the same hot cluster gas creates 
a spectral distortion (see [3] for a review). In high-mass 
clusters (M ~ 10 15 solar masses), the tSZ signal is typi- 
cally a factor of 20 larger than the kSZ signal; however, 
the two signals are comparable for the low-mass clus- 
ters (M ~ 10 13 Af©) which are far more abundant. (For 
brevity, we refer to any object with mass larger than 
10 13 M© as a cluster, even though objects below 10 14 M© 
are usually referred to as “groups”.) The tSZ effect from 
large clusters is now regularly observed in blind surveys 
[4-7], but only upper limits for the kSZ effect from indi- 
vidual galaxy clusters have been achieved to date [8-10]. 

In this paper, we present a clear statistical detection 
of the motions of galaxy clusters through their kSZ sig- 
nal in arcminute-resolution microwave maps made with 
the Atacama Cosinology Telescope (ACT) [11]. Lumi- 
nous galaxies are associated with galaxy clusters [12, 13], 
and we use the Sloan Digital Skv Survey III (SDSS- 
III) Baryon Oscillation Spectroscopic Survey (BOSS) [14] 
catalog of these galaxies as galaxy cluster proxies, giving 
the sky location and redshift for thousands of potential 
clusters. We then treat the microwave temperature mea- 
sured by ACT in the direction of the cluster as a noisy 
estimator of the cluster’s line-of-sight momentum, due 
to the kSZ temperature shift from that cluster. Individ- 
ual cluster momentum measurements have a low signal- 
to-noise ratio, but we combine a large number of dif- 
ferential measurements to obtain estimates of the mean 
relative momentum of cluster pairs in bins of comoving 
cluster separation. This statistic is insensitive to the tSZ 
signal or galaxy emission. The conventional scenario of 
structure formation, driven by gravitational attraction, 
predicts that any pair of clusters should have a slight 
tendency to be moving towards each other rather than 
away from each other [15, 16], and we see the expected 
signal in our data at a statistical significance of 3,8a. 

Surjey Data Sets. We make use of two astronomical 
survey data sets. The first is a 148 GHz sky map from 
ACT, a dedicated microwave survey telescope in the At- 
acama Desert of Chile. The map covers a strip approxi- 
mately 3° wide and 110° long with an angular resolution 


of 1.4', centered on the celestial equator and obtained 
over three observing seasons from 2008 to 2010 [11] . Map 
pixels are 0.5' square, and have a noise per pixel ranging 
from 15 to 25 /iK brightness temperature [17]; the map 
is calibrated to 2% by comparing with WMAP [18]. (A 
similar map at 218 GHz has higher noise and is used for 
Table 1 below.) 

The second data set is a catalog of luminous galax- 
ies from BOSS Data Release 9 (a combination of the 
CM ASS and LOZ samples from DR9), a component of 
the Sloan Digital Sky Survey III [19-21]. The catalog 
contains 27291 galaxies in a 220-square-degree region 
overlapping the ACT sky region (right ascension range 
-43° to +45°). Galaxies are selected to lie at least 1' 
away from any radio source in the 1.4 GHz FIRST ra- 
dio catalog [22]; radio contamination is not a significant 
issue. Spectroscopic redshifts range from z = 0.05 to 
z = 0.8 with a mean redshift of 0.51; luminosities are 
estimated as in Ref. [24] . A halo-model correlation func- 
tion analysis shows that most of the BOSS galaxies re- 
side in haloes with masses around 10 x3 solar masses, with 
around 10% to 15% in haloes as large as 10 14 M© [23]. 

To estimate the microwave temperature distortion T* 
associated with galaxy i , we match-filter the ACT 148 
GHz map with a characteristic filter scale at the map 
resolution of 1.4' [24], to suppress noise from the pri- 
mary microwave background fluctuations. We then fol- 
low the procedure used in Ref. [24] : a 10' by 10' submap 
centered on the galaxy is repixellized into 0.0625' subpix- 
els, convolved with the ACT beam profile to smooth the 
map, and then averaged over all subpixels within 1' of the 
galaxy. The 1' binning radius maximizes the signal-noise 
ratio of our kSZ detection, but varving the binning radius 
between 4" and 4' only changes the detection significance 
by 0.5cr. 

Luminosity correlates with halo mass for the galaxies 
in our catalog. To confirm this, we divide our sample 
into five luminosity bins; Table 1 displays the mean cen- 
tral temperature distortion corresponding to the galaxies 
in each bin. The rightmost column gives the linear com- 
bination of ACT 148 GHz and 218 GHz signals which 
corresponds to the tSZ distortion brightness temperature 
at 148 GHz, with the next-largest component due to dust 
emission projected out [25]. Galaxies in the three highest 
luminosity bins, corresponding to about 20% of the total, 
show mean temperature decrements consistent with halo 
model cluster masses [23] and with the mean temperature 
decrements found in Ref. [24]. 

Mean Pairwise Momentum. Combining the above sur- 
vey data provides a set of galaxy cluster sky positions 
and redshifts (from the luminous galaxy positions and 
redshifts) and line-of-sight momenta (from the ACT tem- 
perature). To compare this data set with cosmological 
models, consider the mean pairwise momentum statistic: 

PpairM — ((Pi “Pi) ' *ij) > 
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TABLE I: Mean brightness temperature fluctuations in sky directions corresponding to the BOSS DR9 galaxies, in bins of 
galax\ luminosity. The right column corresponds to the tSZ brightness temperature at 148 GHz, 5Ttsz = 6 T\ 4 b — 0.325 <57218, 
projecting out a dust emission component [25]. The ACT maps are match-filtered at an angular scale of 1.4', equal to the beam 
size at 148 GHz, then subpixelized, convolved with the beam profile, and summed over all subpixels within 4" of the galaxy. 


where galaxy cluster i has momentum p* and comoving 
position r*, the comoving separation vector between a 
pair of clusters i and j is = r* — r j, overhats denote 
unit vectors, and the average on the right side of the 
equation is over all cluster pairs in a bin around comoving 
separation r = |ry|. If two galaxy clusters are moving 
towards each other, their contribution to PpairM will be 
negative, and if moving apart, positive. An estimator of 
Ppair( r ) using only line-of-sight momenta is [26] 


PpairCr) = 

°ij — r ij • 
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where 6 is the angular separation between two clusters on 
the sky and r* = [r*| is the comoving distance to cluster 
i , which can be computed from the duster redshift us- 
ing standard ACDM cosmological parameters [27]. (The 
cluster velocity gives negligible contribution to the dis- 
tance estimate for clusters at a cosmological distance.) 
The statistic p pa irO") Is equal to the familiar mean pair- 
wise velodty Upair(r) [28-30] times the average mass of 
the clusters in the sample. 

We can measure the line-of-sight component of the mo- 
mentum via the kSZ microwave temperature fluctuation, 
7ksz,i ^ — iVksz assuming that the ratio of the total 

cluster mass to its mass in hot gas is simply the univer- 
sal ratio of matter density to baryon density Q& [31]. 
The normalization Ny^sz depends on the pixel scale and 
beam size of the microwave map, and the cluster density 
profile. Simulations including these effects [32] give an 
expected mean temperature signal in the ACT 148 GHz 
map of 2.2 and 0.9 fiK for clusters with a typical line- 
of-sight velocity of 200 km/s and masses 10 14 and 10 13 
M 0 . 

The statistic Ppair(r) is both linear and differential, 
giving it desirable systematic error properties [33]. Any 
microwave temperature signal associated with individual 
galaxy clusters, like the tSZ effect, will average to zero 


as long as it does not depend on the relative distance 
between cluster pairs. Redshift-dependent signals can 
contribute to p pa ir(r) and be confused with the cluster 
kSZ signal, including infrared emission from galaxies in 
the cluster (which increases with redshift out to z = 2), 
any radio source emission, and small variations of the 
tSZ signal due to evolution of average cluster mass and 
temperature. However, we can measure these, effects on 
average by simply finding the average microwave tem- 
perature T(z) corresponding to clusters at a given red- 
shift, and correcting the temperature in the direction of 
an individual galaxy duster for this redshift-dependent 
piece. We evaluate a smoothed T{z) by averaging the 
temperature towards all galaxies, each with redshift Z\ 
and a gaussian weight factor exp[— (z — z*) 2 /2<t 2 ] with 
<j z — 0.01; our results are nearly insensitive to the value 
of <j z within a wide range. The resulting T (z) has a mean 
near zero and an absolute value of up to 3 fiK. 

We thus evaluate the mean pairwise kSZ signal, cor- 
recting for possible redshift-dependent temperature con- 
tributions, as 


, x m - T(*)) - (Tj - T(Zj )) J Cy _ 

Pksz(r) = ^ — -3 • ( 4 ) 

c ij 


This quantity differs from Eq. (2) by the amplitude fac- 
tor N kSZ . Figure 1 displays this statistic for the ACT 
pixel temperatures corresponding to the 7500 most lumi- 
nous BOSS DR9 galaxies in the ACT sky region (L > 
6.3 x 10 10 L©); this luminosity cut minimizes the total 
noise from combined Poisson and pixel noise. Also dis- 
played is the signal extracted from a kSZ-only sky sim- 
ulation, based on underlying large- volume cosmological 
simulations [32], adjusting the mass limit of the simu- 
lation halos to give the best fit to the data. We infer 
that our galaxy luminosity cut corresponds to a cluster 
halo mass limit of roughly M 200 — 2.5 x 10 13 M© and a 
mean cluster halo mass of M 200 = 4.0 x 10 13 Af©. Er- 
ror bars are estimated via bootstrap resampling. Neigh- 
boring bins have correlated error bars of around 20% as 
determined using independent simulation volumes. 
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FIG. 1: The upper panel shows the mean pairwise momentum 
estimator, Eq. (4), for the 7500 most luminous BOSS DR9 
galaxies within the ACT sky region (red points), with boot- 
strap errors. The solid line is derived from numerical kSZ sim- 
ulations [32] using a halo mass cutoff of M 200 = 2.5 X 10 13 M©. 
The measured signal is 3.8<r away from zero including bin cor- 
relations. The lower panel displays the same sum but with 
randomized map positions, and is consistent with null signal, 
0.4(7 away from zero. 


The measured points largely fall below zero and are 
consistent with the simulated signal; they differ from a 
null signal by 3.8<r including correlations. The measured 
points approach zero signal as the comoving pair sepa- 
ration increases, which demonstrates that the signal de- 
pends on spatial separation, not redshift separation. 

Null tests are simple, as the statistic is essentially a 
sum of pixel temperatures, half with positive and half 
with negative signs, with weights corresponding to rel- 
ative galaxy positions. Figure 1 also displays the null 
test corresponding to using the same weights but random 
pair positions compared to the signal plot (Ax 2 = 11.6 
for 15 degrees of freedom). Success of this null test veri- 
fies that the function T ( z ) correctly models any redshift- 
dependent contributions to the microwave signal. Chang- 
ing the sign in the second term of Eq. (4) from negative to 
positive also gives a null signal (Ax 2 = 9.9 for 15 degrees 
of freedom). 

Discussion and Prospects. The signal in Fig. 1 repre- 
sents the first measurement of the cosmic velocity field 
made directly with respect to the rest frame of the uni- 
verse. It is consistent with simulations based on the stan- 
dard cosmological model. This signal is also the first 
detection of the kinematic Sunyaev-Zel’dovich effect. A 


recent attempt by Kashlinsky et al. to measure the large- 
scale bulk flow via the galaxy cluster kSZ signal uses 
galaxy clusters from X-ray surveys and searches for an 
overall dipole dependence of the microwave temperature 
in the WMAP data at these locations [34, 35]. How- 
ever, Keisler [36] found the first reported detection was 
not statistically significant. Osborne et al. [37] reana- 
lyzed the most recent results including both a monopole 
and dipole term, obtaining limits on a bulk flow a factor 
of three below the reported detection of Ref. [35] . Mody 
and Hajian [38] also fail to reproduce the bulk flow result 
using Planck and ROSAT galaxy clusters. 

The detection of a non-zero mean pairwise momen- 
tum from a kSZ signal presented here can also be inter- 
preted as a measure of baryons on cluster length scales; 
a deficit of observed baryons has long been a cosmologi- 
cal puzzle [39]. Our signal is roughly consistent with the 
standard baryon fraction based on primordial nucleosyn- 
thesis, given independent halo mass estimates based on 
clustering of our luminous galaxy sample. This issue will 
be addressed in a future paper. 

Future improved measurements of the mean pairwise 
velocity have the potential to put strong constraints on 
dark energy and modified gravity [40-42] . The measure- 
ment we have presented here is the first step on a new 
path to constraining structure growth in the universe. 
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